Bose-Einstein Condensation of Metastable Helium 
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We have observed a Bose-Einstein condensate in a dilute gas of 4 He in the 2 3 Si metastable state. 
We find a critical temperature of (4.7 ± 0.5) /iK and a typical number of atoms at the threshold 
of 8 x 10 6 . The maximum number of atoms in our condensate is about 5 x 10 5 . An approximate 
value for the scattering length a = (16 =b 8) nm is measured. The mean elastic collision rate is then 
estimated to be about 2 x 10 4 s _1 , indicating that we are deeply in the hydrodynamic regime. The 
typical decay time of the condensate is 2 s, which places an upper bound on the rate constants for 
2-body and 3-body inelastic collisions. 
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Bose-Einstein condensation (BEC) of dilute atomic 
gases was first observed in alkali atoms in 1995 and then, 
a few years later in atomic hydrogen. Since then, the field 
has developed in a spectacular way both experimentally 
and theoretically |]] . So far only condensates with atoms 
in their electronic ground state have been produced. 

Several laboratories are currently involved in the 
search for BEC of atoms in an excited state, namely no- 
ble gases in an excited metastable state. Helium in its 
triplet metastable 2 3 5i state ( 4 He*) is of particular in- 
terest. The first advantage of 4 He* is its large internal 
energy (19.8 eV). It allows for a very efficient detection of 
the atom by ionization after collision with another atom 
or a surface, which can be of interest for atomic lithog- 
raphy [§J3). Second, helium is a relatively simple atom 
which allows for quasi-exact calculations that are useful 
in metrological applications. Third, mixtures of 3 He and 
4 He can be used to study quantum degenerate mixtures 
of bosons and fermions. Finally, Penning collisions are 
expected to be inhibited for spin polarized atoms due to 
spin selection rules. This effect, first pointed out in [[|, 
was confirmed by subsequent calculations ||. 

The present article describes the observation of BEC 
of 4 He* atoms. Similar results have also been obtained 
at IOTA, Orsay Q. The two experiments differ by their 
detection methods. The Orsay group detects the atoms 
falling on a microchannel plate, whereas we use an opti- 
cal absorption imaging of the atomic cloud on a CCD 
camera. The two experiments therefore give different 
and complementary information on the physics of BEC 
in 4 He*. 

The first step of our experiment is the efficient load- 
ing of a magneto-optical trap (MOT). The experimental 
set-up is described in detail in (?]]. A discharge atomic 
source ensures a high flux of triplet metastable atoms of 
10 14 atoms/s.sr, with a mean velocity of about 1000 m/s. 
The atomic beam is collimated || and Zeeman slowed by 
laser light at 1083 nm (2 3 5i-2 3 P2 transition). A narrow 



frequency band master oscillator (DBR diode laser) in- 
jects a Yb-doped fiber amplifier with an output power of 
500 mW. Using this set-up, it is possible to trap ~ 8 x 10 8 
atoms in the MOT at a temperature of 1 mK. 
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FIG. 1. Top view of the magnetic trap. Coils Qi and Q2 pro- 
duce a quadrupole field used for the MOT. Combined with a third 
coil C3, they produce a magnetic IofTe- Pritchard trap. Helmholtz 
coils Hi and H2 compensate the bias field. The probe beam used 
for optical detection is perpendicular to the longitudinal axis of the 
trap. 

4 He* atoms are confined at the center of a small 
(4x4x5 cm) quartz cell. All coils are external to the 
cell (see Fig.|l|). The coils Qi and Q 2 (144 turns and 7 
cm diameter) combined with the coil C3 (108 turns and 
4 cm diameter) produce an anisotropic magnetic Ioffe- 
Pr it chard trap. Two additional Helmholtz coils reduce 
the bias field, in order to increase the radial confinement 
of the trap. A current of 45 A in all the coils produces 
a 4.2 G bias field, radial gradients of 280 G/cm and an 
axial curvature of 200 G/cm 2 . These values correspond 
to trapping frequencies of 115 Hz in the axial and 1090 
Hz in the radial directions. The current in the coils can 
be switched off in 200 /is. 

The second step of the experiment is the loading of 
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the magnetic trap. After switching off the MOT field, 
the cloud is further cooled down to about 300 jaK during 
a 1 ms optical molasses phase. To increase the trans- 
fer efficiency from the molasses to the magnetic trap, 
the atoms are optically pumped by a circularly polarized 
laser pulse. 3 x 10 8 atoms are loaded in the magnetic 
trap. The lifetime of the atomic cloud in the magnetic 
trap is about 35 seconds, and its temperature is 1.2 mK 
after compression and bias compensation. 

The last step of the experiment consists of evaporative 
cooling performed by radio frequency (RF) induced spin 
flips. The frequency is ramped down from 160 MHz to 
around 12 MHz in 15 seconds. After evaporation, the 
trap is switched off and the cloud released from the trap 
is probed by absorption imaging on a CCD camera whose 
quantum efficiency is 1.5% at 1083 nm. 



The spatial distribution of the absorption pictures is 
fitted with the sum of two functions, one for the conden- 
sate and one for the thermal cloud (see insert in fig. ||). 
The function for the condensate is an integration along 
the z-axis of a paraboloidal distribution. It describes 
the equilibrium density profile of the condensate within 
the harmonic trap in the Thomas- Fermi limit ]l0| . The 
function for the thermal cloud is a function valid for 
a bosonic gas close to the transition where the chemi- 
cal potential fi is an adjustable parameter [Q. From the 
fit, we extract the ratio between the number of atoms in 
the condensate, TVo, and the total number of atoms, N. 
Plotting No/N versus the temperature T (Fig.||) gives 
the value of the critical temperature T c = 4.7 ± 0.5 /iK, 
which is confirmed by a TOF measurement performed on 
a thermal cloud just above the transition. 
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FIG. 2. a) Time-of-flight absorption images (1.4 x 1.4 mm) after 
an expansion time of 0.2, 2 and 8 ms; b) Ellipticity of the conden- 
sate for increasing expansion times. The solid line is the theoret- 
ical prediction without any adjustable parameters. The inversion 
of ellipticity is characteristic of the behaviour of an expanding con- 
densate. 

When the RF frequency of the evaporation is ramped 
down to a final frequency below 13 MHz, a narrow struc- 
ture appears on the absorption image which we identify 
as a condensate. The strongest evidence for the pres- 
ence of the condensate is the evolution of the shape of 
this structure when released from the trap in the time- 
of-flight (TOF) measurement. Its anisotropy increases 
as the expansion time is increased, and its ellipticity un- 
dergoes an inversion (see Fig.|2| a). This observation is a 
consequence of the mean-field interaction between atoms 
in the condensate ||. The theoretical prediction con- 
taining no adjustable parameters agrees well with our 
measurements (see Fig.^Jb). 
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FIG. 3. Condensed fraction No/N versus T. The insert shows a 
ID profile of an absorption image, displaying a bimodal structure, 
composed of a condensate and a thermal component, which is char- 
acteristic of the bosonic gas below T c . The thermal component was 
fitted for clarity. 

The results presented until now are based only on mea- 
sured sizes of the clouds and on relative numbers like 
No/N which do not require absolute calibration of the 
optical detection. The procedure described previously 
for the measurement of ellipticity and critical tempera- 
ture is not suitable for absolute measures of N and Nq. 
So, we use a different procedure to calibrate the num- 
ber of atoms at threshold. The current is switched off 
in two steps: first, in the Helmholtz coils, then, after a 
delay of 10 ms, in the 3 coils Qi, Q2 and C3, whereas 
the current in all five coils was switched off simultane- 
ously for the measurements above. Then, we take an 
absorption image with an exposure time T exp = 200 /is 
after a TOF time of 5 ms. This procedure minimizes 
the effects of eddy currents of the Helmholtz coils which 
last several milliseconds, Zeeman shifting the atomic res- 
onance with respect to the probe beam frequency Q. 
There are two other obvious sources of sensitivity losses. 
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First, the width of the absorption lineshape is measured 
to be twice the natural linewidth, which we attribute 
to laser linewidth and power broadening. This reduces 
the absorption cross- section by a factor 2. Second, the 
probe beam propagates perpendicularly to the x-axis of 
the Ioffe-Pritchard trap (see Fig.|l|) along which atoms 
are polarized. It would therefore be necessary to take 
into account the populations of the different Zeeman sub- 
levels, and the Clebsch-Gordan coefficients of the various 
Zeeman optical transitions excited by the probe beam. 
But any residual magnetic field redistributes the popu- 
lations among Zeeman sublevels. Assuming equally pop- 
ulated sublevels gives an extra loss by a factor 9/5. Fi- 
nally, there is another source of losses which is specific to 
4 He* when optically detected and particularly important 
at the high densities obtained in our experiment. The 
metastable atoms have a huge Penning ionization cross- 
section in the presence of resonant light, so that losses 
can accumulate during the probe pulse at large densi- 
ties. Unfortunately, probing at low intensity is no longer 
possible if T exp is made shorter because of the low ef- 
ficiency of the CCD camera. This long exposure time 
also increases the acceleration of the atoms due to radi- 
ation pressure which pushes them out of resonance [ [l2| . 
In the present stage of the experiment, it is difficult to 
give a quantitative description of the combined effect of 
these phenomena. So, we prefer to wait an expansion 
time of 5 ms to sufficiently reduce the atomic density. 
Indeed, we observe that the measured number of atoms 
is an increasing function of the expansion time reaching 
a plateau after 5-6 ms. After the corrections mentioned 
above, the total number of atoms N c at the threshold is 
measured to be 5xl0 6 atoms with an accuracy of about 
50%. 

N c can also be estimated from N c = 1. 202 (k bT c / 'ho) 3 , 
where Co is the geometrical average of the frequencies of 
the trap Jl3|. We deduce N c = 8.2 x 10 6 atoms with an 
uncertainty of about 30% compatible with the previous 
value. Because it is easier to trace down the error on 
N c when derived from T c , we arbitrarly choose this value 
of N c to estimate the scattering length a. Also, assum- 
ing that the transition occurs at a phase space density 
equal to n(0)A^ B = 2.612, where Xjb is the de Broglie 
wavelength, we can also derive the density at the center 
of the trap n(0) = (3.8 ± 0.7) x 10 13 atoms/cm 3 at the 
transition. 

Within the Thomas-Fermi approximation, one can ex- 
tract the chemical potential \i from the size of the con- 
densate |l(J. As the optical detection around the tran- 
sition has been calibrated, we can now use our previous 
relative measurements to deduce the absolute number of 
condensed atoms No below the transition. Typical values 
of /i = 1.4 x 10~ 29 J and N = (4 ±1.5) x 10 5 are obtained 
with condensates prepared at temperatures ranging from 
1.2 /iK to 3 jaK. Finally an estimation of the scattering 
length a can be given using a = a/lbNo x (2/i/?kZ>) 5 / 2 , 



where a = (Ti/muu) 1 / 2 is the characteristic size of the 
ground state of the trap. We find a = (16 ± 8) nm which 
is compatible with the value given by recent theoretical 
works |||. The error is mostly due to the uncertainty 
on the number of atoms. Knowing the scattering length 
a, the density at the center of the trap n(0) and the 
critical temperature T c , we obtain a mean rate of elas- 
tic collisions % Q u ~ 2 x 10 4 s _1 near threshold, leading 
to Cu/jcoii = 0.17. We thus enter in the hydrodynamic 
regime, an interesting feature for a gas above T c rB~S|- 
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FIG. 4. Decay of the condensate. This measurement was per- 
formed with an RF shield at 12.3 MHz. 

Fig.|| shows the evolution of the number of atoms in 
the condensate versus the trapping time. The typical 
lifetime is about 2 s. In the present stage of the ex- 
periment, it is not possible to discriminate between 2- 
body or 3-body decay. However, assuming that only 2- 
body collisions lead to losses in the condensate, we can 
place an upper bound on the collision rate constant G 
between spin polarized 4 He* atoms, which is predicted 
to be inhibited by a factor 10 4 compared with the rate 
constant for Penning ionization collisions between unpo- 
larized atoms. The evolution equation of the number of 
atoms N(t) = — G J n 2 (r)d 3 r, integrated on the whole 
condensate in the Thomas-Fermi approximation, gives 
the evolution of the number of atoms Nq in the conden- 
sate. A fit leads to G < (4.2 ± 0.6) x 10" 14 cm 3 /s, which 
corresponds to a reduction factor larger than 2 x 10 3 , 
much larger than the previously measured ones ]l9| , p0[ , 
in agreement with theoretical calculations [QJ^]. 

If one assumes now that 3-body collisions are respon- 
sible for the decay of TVo, one can as well give an up- 
per bound for the rate constant L defined by N(t) = 
-L J n 3 (r)d 3 r. Fitting our data gives L < (2.8 ± 0.2) x 
10 -27 cm 6 /s. This value can be compared to theoretical 
predictions pl|,22[ : for example, our upper limit is com- 
patible with |2l[, which finds L = 3.9fta 4 /2ra = 2 x 10~ 27 
cm 6 /s. 

In conclusion, this article shows the evidence for the 
formation of a BEC of helium atoms in the metastable 



3 



state 2 3 5i. Our results concerning the losses due 
to inelastic collisions show that the spin polarization 
does inhibit the Penning ionization collisions between 2 
metastable helium atoms by more than 3 orders of mag- 
nitude, as theoretically predicted. We find a large scat- 
tering length, which results in very large rates of elastic 
collisions : the cold gas at threshold is in the hydrody- 
namic regime, which we plan to study in more details. 
We also plan to further exploit the original characteris- 
tics of this new born condensate of atoms in an excited 
state. 
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